A number filter for matter-waves 
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In current Bose-Einstein condensate experiments, the shot-to-shot variation of atom number fluctuates up 
to 10%. In here, we present a procedure to suppress such fluctuations by using a nonlinear p—iz — p matter 
wave interferometer for a Bose-Einstein condensate with two internal states and a high beam-splitter asymmetry 
(p,p 7^ 0.5). We analyze the situation for an inhomogeneous trap within the Gross-Pitaevskii mean-field theory, 
as well as a quantum mechanical Josephson model, which addresses complementary aspects of the problem and 
agrees well otherwise. 
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I. INTRODUCTION 

Nonlinear optical wave propagation is known to give rise 
to chaos, spectral and temporal distortion or an amplifica- 
tion of noise This fundamentally limits the signal- 
to-noise ratio in state-preparation experiments or measure- 
ments. For example, the nonlinearity constrains the capac- 
ity of optical communication systems |2], or the resolution 
of a gravitational- wave interferometer [3|], where the momen- 
tum transfer to the mirror produces an intensity-dependent 
phase shift. However, optical nonlinearities are also capable 
of wave-packet self-stabilization and of a phase-sensitive re- 
duction of noise. Second- and third-order nonlinearities, and 
especially the Kerr effect in optical fibers, produce energy sta- 
bilization and a noise reduction below the standard quantum 
limit iHJlSlllllllIllIl in various experimental configu- 
rations IU2I1 . 

The past decade of matter-wave physics has also shown re- 
markable similarities with the development of quantum optics 
in the 60's. A lucent description of this parallelism of quantum 
optics [13] and atomic matter waves is found in 11411 . Starting 
from the seminal measurement of spatial coherence in normal 
and degenerate gases 11811 . the field has eventu- 

ally progressed to study density fluctuations in trapped, three 
dimensional Bose-Einstein condensates (BEC) 11191 l20l 12111 
and fermionic gases [22]. By reducing dimensionality via ge- 
ometric confinement in planar traps, one-dimensional traps, 
in optical lattices or on atomic chips Il23ll the field has now 
been opened to a plethora of condensed matter phenomena 
SSESIlSISIlSIMIll. While the perfection of 
communication quality is the key issue for optics today, the 
main application for cold atomic matter waves is quantum 
metrology and sensing. Reaching the quantum limit and sur- 
passing it with matter-waves is a major research direction 
lHHlHmmiiSmiiiEIOIii. In this con- 
text, atoms or ions prove to be more flexible than light, as we 
can control many-particle entanglement and exploit different 
quantum statistics JH H3, EH . 

The statistical ensembles that are generated in most of the 
current experiments are never of the quality as theoretically 
envisaged. In particular, most of the current BEC experi- 



ments face a shot-to-shot variation of particle number N of 
about 10%. This is primarily due to technical uncertainties in 
the evaporation procedure. If each individual BEC realization 
would be characterized by a pure Fock state I^Pw), then the 
uncertainty &n m atom number will lead to an mixed state 
ensemble with a density operator 



p = £^|¥ w ) CIV | 

N 



(1) 



Thus, each observable will lose contrast caused by this num- 
ber uncertainty. In this paper, we will establish an atom num- 
ber filter for matter waves that allows a number stabilization, 
i. e., after passing each BEC through the filter the number un- 
certainty is less than before. 
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FIG. 1: Setup for an asymmetric optical nonlinear interferometer 
with a propagation length L and a splitting ratio p : ( 1 — p), p ^= 0.5. 
Due to a Kerr-nonlinearity (susceptibility a), one obtains a differ- 
ential phase shift 0nl = 01 — <h = a i^-P — l)Am> proportional to an 
input photon intensity N[ a . A subsequent self-interference of strong 
a i and weak field «2 stabilizes the output field intensity N ou t = \ a out | 2 
and diverts the noise to the rejection port a le j. 



This can be achieved by using a nonlinear matter wave in- 
terferometer, cf. Figs. [T]and [3] which is in analogy to a nonlin- 
ear fiber optics setup [7[. We will assume that the condensate 
consists of atoms with two internal states. Due to the highly 
asymmetric splitting, which is crucial in this setup, the con- 
densate fraction in one arm of the interferometer experiences 
a strong nonlinear phase evolution, while the other part only 
undergoes a weak nonlinear phase shift. 
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The underlying physical mechanism of the suppression 
of number fluctuations is based on the repulsive interaction 
amongst particles and has been used in the context of spin 
squeezing @|47tl or the Josephson effect H [H, [U . The 
ideas presented in here are also related to the work of Poulsen 
and M0lmer lf5lll . since their approach combines ideas for 
light squeezing in a nonlinear optical interferometer with the 
idea for spin-squeezing of two spatial, initially identically oc- 
cupied, condensate modes, generated via Bragg scattering. 
Our approach is different, as it explicitly requires a highly 
asymmetric splitting p ^ 0.5 or would it disappear at all and 
it uses internal states of the atom. 

This paper is organized as follows: First, we briefly review 
the central idea of amplitude stabilization of a nonlinear opti- 
cal interferometer in Sec. [II] Second, we introduce an equiva- 
lent model for a bosonic matter wave in Sec. [Ill] This is stud- 
ied within a mean-field picture to consider the effects of inho- 
mogeneous traps as well as a Josephson model of two quan- 
tized plane wave modes, which address the quantum aspects 
and effects of finite particle numbers. Finally, conclusions are 
drawn in Sec. [TV] 



H. THE PRINCIPLE OF NONLINEAR 
INTERFEROMETERS IN OPTICAL FIBERS 

A very fundamental type of nonlinearity, which is present in 
many systems, is the intensity-dependent phase shift. In pho- 
ton optics it is due to the optical Kerr effect [1] characterized 
by a susceptibility a and in matter waves it is caused by inter- 
atomic atom forces measured by the s-wave scattering length 
a s . This leads to a self- or cross-phase modulation and pos- 
sibly to a self-trapping potential in the nonlinear Schrodinger 
equation of motion of wave-packets. 

In the context of nonlinear interferometry Q7|] , the intensity 
filtering property is best in a highly asymmetric, highly trans- 
missive configuration, depicted in Fig.Q] The interference of 
a strong wave with a weak wave can eliminate a major fraction 
of the input noise of a„, in the output port a out . A predomi- 
nant part of the noise is channeled to the rejection port a re j, 
consuming a small fraction of the input photon number. 

Its two-step nonlinear self-stabilization mechanism is very 
simple and visualized in Fig. [2] After the first beam-splitter, 
the asymmetric splitting of an input beam causes an intensity- 
dependent differential phase shift between the two arms a\ 
and «2 of the interferometer. The nonlinear phase shift trans- 
forms intensity amplitude increase/reduction due to field fluc- 
tuations into a correlated phase advance/delay. Therefore, it 
causes a correlated phase spread relative to the average non- 
linear phase shift. The second beam-splitter superposes both 
interferometer beams coherently. It changes the quadrature 
angle of the field amplitude relative to the noise distribution so 
that the amplitude-phase correlation eliminates the amplitude 
noise to a large degree. For perfect stabilization, the phase 
advance/delay of the stronger mode relative to the weaker, 
quasi stationary, linearly propagating mode reduces/increases 
the output transmission by the right amount to eliminate the 
input intensity fluctuations. 



It has been shown that the stabilization mechanism does not 
only eliminate classical noise or works only with continuous- 
wave coherent light. Instead, this method is also applicable 
with broadband ultrashort solitons and in the quantum regime 
of field fluctuations. The fiber-optic asymmetric Sagnac inter- 
ferometer has been used as a photon number filter for optical 
solitons 

QUI- 

Some of the best squeezing results have been 
obtained with this set-up, which did not require any active 
stabilization. The quantum noise reduction below the shot- 
noise (Poisson) limit has been modeled by the quantum non- 
linear Schrodinger equation (NLSE). It is in perfect agreement 
within the measurement uncertainty and stability has been ob- 
tained. Again, the noise reduction mechanism can be readily 
understood by modeling the essentials of the asymmetric in- 
terferometer by linearized fluctuations in a semiclassical ap- 
proach where now the uncorrected vacuum fluctuations are 
entering through the unused input port of the interferometer 
and where the associated phase is the soliton envelope phase. 
The corresponding semiclassical picture is then well repre- 
sented by Fig. [2] where the noise ellipses are then the mini- 
mum uncertainty regions. 




FIG. 2: Schematic representation of the amplitude stabilization 
mechanism of a m in an asymmetric, nonlinear interferometer in 
phase space with quadratures components (a c ,a s ). After the first 
beam-splitter (1), the stronger field a\ experiences a large nonlin- 
ear phase shift that translates amplitude fluctuations into correlated 
phase delay. One can neglect the nonlinear phase shift of the weaker 
field ai- In a second beam-splitter (2), a% is coherently added to the 
phase-shifted field a\ to cause output in a out . This stabilizes the out- 
put fluctuations well below the input level Aa; n 2> Aa out , as indicated 
by the noise ellipses. 

Thus, the question arises, how the analogy of interference 
of bosonic fields can be applied to matter waves. The analogy 
is not trivial, because the quanta of the optical field and the 
massive bosons of the matter field are described by different 
ensembles. Also, on a practical side, it can be asked how well 
the asymmetric interferometer can function as a number filter 
for matter waves. 



IH. MODELING A NONLINEAR INTERFEROMETER 
WITH BOSONIC MATTER- WAVES 

Let us consider a trapped BEC consisting of two-level 
atoms labeled by a — e,g. The complete quantum states are 
then denoted by |c7,k,j) with the internal state a and momen- 
tum component k<j. The possibly time-dependent trapping 
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potentials for the two species are V a (r,f). They are identical 
V e (r,t) — V g (r,t) +A up to a detuning A. The two states are 
coupled by a classical traveling laser field £2(f )<?' kr with the 
time-dependent Rabi frequency Q.(t ) and the wave vector k. 
This configuration represents a Ramsey-Borde-interferometer 
(see Fig. [3} in the standard setup of atom interferometry f52ll . 
Initially, the BEC is prepared in the \g) -state and at an instant 
t = a p-pulse creates a superposition of \g) and \e) with a 
splitting ratio of the populations of p : (1 — p). After a time 
delay T, a further ft -pulse gives rise to an inversion of the pop- 
ulations. After another time interval T, a final beam-splitter 
with splitting ratio p : (1 — p) mixes the populations again. 



31 P \e,k+k) 

I ft K> 




FIG. 3: Setup for a matter-wave interferometer. Absorption of a pho- 
ton k implements an asymmetric beam-splitter with splitting ratio 
p : (1 — p). After a free time evolution of duration T, an optional 
7r-pulse inverts the populations. The second beam-splitter with split- 
ting ratio p : (1 — p) mixes the states for the final detection of one 
channel and the rejection of the other (comp. Fig.QJ. 



A. Theoretical description 

In principle, the dynamics of this system is governed by the 
Schrodinger equation ihd t |y/(f)) = "(0 I V(0) f° r a man y- 
particle state |y/(f)) m Fock space and with the Hamiltonian 



where 



o=e,g 



n 2 v 2 



(2) 



2m 



+ V a (r,t) 



d a (r), 



V d (t) = / d 3 r [ftn(f)e ikr 4(r)<^(r) +h.c] , 



+8ggal(r)a f g (r)a g (r)a g (r) 
+2g eg dl(r)al{r)d g {r)a e (r)] . 

In there, we introduced bosonic field operators 
[a M (x),4(y)] = <5(x-y)5 MV 



(3) 



and interatomic coupling constants g^ v = AntP'a^/m be- 
tween same and different species. a^ v is the corresponding 
scattering length and m is the mass of a single particle. 



In order to get physical insight into this problem, we con- 
sider two simplified scenarios. On the one hand, we study the 
classical field approximation in Sec. lIIIBl where operators are 
replaced by complex amplitudes a a — * a a , which yields the 
Gross-Pitaevskii (GP) equation. Thus, the statistical character 
of the operator is neglected. Furthermore, we will consider a 
quasi one-dimensional, cigar-shaped configuration with tight 
confinement in the radial direction [30]. The radial part can 
then be integrated out directly, which results in modified cou- 
pling constants g^ v — * g^ v . In the following, we tacitly drop 
the bar. On the other hand, the operator character is accounted 
for in the Josephson approximation of Sec. lIIICl In there, we 
neglect spatial inhomogeneity and consider only the behavior 
of two plane wave modes. 



B. Classical field approximation 

Within the classical field approximation, the corresponding 
scaled quasi-one-dimensional GP equation reads 



a g {z,t) 



Hgp(0 = -R^ + 



a,. 



(4) 



V e (z,t)+v e n(f)e' te 
Q.*(t)e-^ V g (z,t) + v g 



where the mean-field energies are v M = g^n^ +^v«v with 
ju 7^ v S {e,g} andncj = \a a {z,t)\ 2 . 

In the general case of time-dependent pulses and spatially 
inhomogeneous traps, it is only possible to solve this equa- 
tion numerically. Results of such a calculation are presented 
in Sec. lIIIB"2l However, if the system size is large and time- 
dependent pulses happen on short times, then one can solve 
this simplified situation analytically and gain qualitative un- 
derstanding. 



1. A bulk BEC in the Raman-Nath approximation 

The beam-splitter is realized by a quasi-instantaneous (t <C 
T, 1/A) p-pulse in the form of a traveling laser wave with a 
Rabi frequency Q.. Mathematically, this is described via the 
unitary transformation 



Up = e -ie(e^ ++ Kc.) = 



cos 9 
-/sin 9e~' kz 



/sin 0e' kz 
cos0 

£2t 



(5) 



In the following, we have propagated the mean-field state 
with flat potentials V g (z,t) — 0, V e (z,t) = A for a total time 
2T in Eq. (0). The free nonlinear evolution was interrupted 
by the interferometer pulse sequence p — % — p depicted in 
Fig. [3] We assume that all the population n is initially in the 
ground state component {a e (0),a g (0)} = {0, ^fne lk z z }, mov- 
ing with a generic momentum k g . The general solution for the 
free propagation is obtained easily by making the plane-wave 
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ansatz 

{a e (z,t),a g (z,t)} = {ccee^-^^e^-^}, (6) 

with time-dependent phases (j) a (t) and a recoil-shifted mo- 
mentum k e = kg + k. We are interested in the particle density 
of the output channel n° ut = n a (2T) = \ a a {2T)\ 2 of the inter- 
ferometer. After simple algebra, one finds for the transmitted 
channels 



n° ut = n(t;-y cos [2Tn82(p-l)]) 



(7) 



$=PP+(i -/>)(! -P), Y=2y/pp(l-p)(l-p), 

with &2 = gee — 2geg + g g g a central difference of scattering 
lengths. The population in the other channel «° ut = n — n° ut 
follows from number conservation. 

A simplified but less efficient form of the nonlinear interfer- 
ometer is found from a symmetric mixing p = i at the final 
output beam-splitter 



- V/X 1 -p)cos \2Tn&i{p - \)] 



(8) 



The most salient features of Eqs. (O and © are the absence 
of linear phase shifts due to the intermediate 7T-pulse and the 
nonlinear phase shift. It is proportional to the total interaction 
time 2T, to the density n, the central difference 52 of scatter- 
ing lengths, as well as the asymmetry p ^ 0.5 of initial beam- 
splitting. The oscillatory response of the interferometer with 
respect to a varying input particle number n stabilizes the out- 
put particle number n° ut , if operated in the vicinity of an ex- 
ternum. This suppression of number fluctuations represents 
a nonlinear number filter for matter waves and is depicted in 
Fig. 2] It is also important to note that a symmetric splitting 
p = j, or vanishing central difference 52 = of scattering 
lengths lead to no effect at all. 

Alternatively, if one considers a simplified interferometer 
setup in Fig. [5J i. e., without the intermediate 7T-pulse and 
chooses a symmetric final beam-splitting p = \, one obtains 



= n ( \ + ^p(l-p)cos[{A D + (5i + p&i)n}T] 



(9) 



with a total propagation time T. This interferometer is sen- 
sitive to single particle phase shifts, here in particular to the 
Doppler-shifted detuning Ad = A + k\ /2 and also the differ- 
ence of scattering lengths 5i = g eg — g PP . For example, this 
setup is used for atom gravitometry 15211 . but it is also more 
susceptible to experimental noise (frequency jigger), which 
deteriorates visibility. Nevertheless, it is favorable for 87 Rb 
BEC's, where the coupling constants are g ee '■ geg '■ ggg = 103 : 
1 : 0.97, hence 52 vanishes. In contrast, there is an effect in 
the latter setup, since the difference 5i is finite. 



2. 



An inhomogeneous BEC in a square well trap 



In order to study the reduction of the filter performance 
caused by the inherent inhomogeneity of trapped atomic 
BECs, we have chosen a square-well potential 



V g (z,t) 



- ) v s 



Vg < 0, t < 0, \z\ <L 



0. 



t >0 



(10) 



to account for the initial inhomogeneity of the ground-state. 
After the first /5-pulse, the trap is switched off permanently. 
The excited state potential is identical, but shifted in energy 
by the detuning from the laser, i. e., V e (z,t) — V g (z,t) + A. 
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FIG. 4: Response of a highly asymmetric p = p = 0.9, nonlinear 
p — n — p matter- wave interferometer, a) Outgoing particles in the e- 
channel N° ut versus the incoming number N' g n in the g-channel: ho- 
mogeneous mean-field (dashed-dotted), inhomogeneous mean-field 
(dotted), and two-mode approximation (solid). The line marked with 
□ is the trivial response of a — 71 — interferometer, i. e., JV° ut = 
and can be used to assess the loss of particles in the unobserved 
channel, b) Corresponding normalized number fluctuations s° ut in 
a semilogarithmic representation. The optimal working point for the 
interferometer as a number-filter device is 2V? = 940 (see Fig. [5j- 
There, one finds a strong, sub-shot noise suppression of number fluc- 
tuations in the output channel s° ut <C 1 . 



In the numerical simulations for the homogeneous and in- 
homogeneous condensates, we have used generic parameters 
for a quasi- Id elongated 87 Rb BEC with an atomic mass 
m = 87amu. In a prolate harmonic oscillator with a trap- 
ping fre quency co z — 4 1/s, the basic length unit would be 
flho = yhfnW) z = 13.2 /im. We have used a multiple of this 
scale for the length of the square-well trap L = 132/xm(10). 
Potential depths V g = — /;60Hz(— 15) and detunings A = 
are measured in natural energy units hco z . In a copropagating 
Raman laser configuration, one can have a vanishing momen- 
tum transfer k = and we assumed that the condensate is ini- 
tially at rest k g — 0. Short laser pulses were used such that the 
beam-splitters were highly asymmetric p = p = 0.9 and the 
propagation time between the pulses was T = 20 ms. We have 
deliberately modified the s-wave scattering lengths for 87 Rb 
slightly to obtain dimensionless quasi- Id coupling constants 
gee = 0.034, geg — 0.10, and g gg = 0.068. Thus, the superior 
p — n — p interferometer scheme remains applicable. In prin- 
ciple, this can be achieved via Feshbach resonances or using 
other elements like 85 Rb or 23 Na altogether. 

For this situation, we have numerically solved the time- 
dependent inhomogeneous GP Eq. and find a similar be- 
havior as for the homogeneous limit in Fig. |4] As expected, 
we obtain a number stabilization of the nonlinear filter, but 
at a slightly reduced performance due to the inhomogeneous 
averaging. We have also verified that the best number stabi- 
lization is achieved for highly asymmetric splitting, that the 
effect does not occur for equal scattering lengths, or linear 
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matter-wave interferometers at all. Moreover, the assumed 
square-well trap is not a peculiarity in this context, as we have 
tried a harmonic oscillator trap and find qualitatively similar 
results. 



C. Quantum mechanical two-mode approximation 

In the classical field approximation for a homogeneous bulk 
system of Sec. IIII B 1 1 we have examined the static number 
filter response of the interferometer. The macroscopically oc- 
cupied amplitudes were described like two coherently cou- 
pled nonlinear oscillators. They exhibited a noise suppression 
that is analogously used in many other physical systems rang- 
ing from electrical circuits to coupled Josephson junctions 

In order to probe the quantum aspects of such a system, we 
will assume now that the atomic fields a a (z) are dominated by 
two plane-wave modes labeled with bosonic field amplitudes 
e and g 



a e {z) 



JkeZ 



a g {z) 



Jk s z 



(11) 



Their residual coupling to other modes 8d a , is small at the 
relevant time scales and will be disregarded altogether. 

Within this Josephson approximation and in the quasi- Id 
configuration, we can simplify the Hamiltonian of Eq. (O fur- 
ther to H w Hj with 



k 2 

Hj = A D n e + -j-flg + £l{t)e'g + Q* (t)g T e+ 
\g ee h e (h e - 1) + \g gg n g (n g - \)+g eg n e n g , 



(12) 



where n a — {e*e ,g T g } denotes the particle number operator 
for the two modes and Ad = A + k 2 /2 is the Doppler-shifted 
detuning as in Eq. (O. The consistency of the limit can be 
checked quickly by replacing the mode operators again by 
complex numbers, just to recover Eq. (|4j. 

Clearly, the Josephson Hamiltonian conserves the particle 
number, 



N = n e + n Hl [Hj,N]=0. 



(13) 



Thus, a general state in the Af -particle sector of Fock space is 
given by a superposition of the states |n e ,n g ) with N = n e +n g 



\y N (t)) = 1 £y%(t)\N-n,n) 

n=Q 



(14) 



The time evolution of such a state leads to a simple one- 
dimensional difference equation for the time-dependent am- 
plitudes 



M(t) =^^ + q"{t)r N - X +q n+l *{t)r+\ (15) 
k 2 

w" =A D (N -n) + j-n + \g gg n{n - 1)+ (16) 
\g ee (N-n){N-n-\)+g eg n{N-n), 



q n {t) =Q(f)V*(#-»+l)> 



(17) 



which can be solved easily on a computer or approximated 
analytically 115611 . 

We are now in a position to discuss the number stabi- 
lization scheme on a full quantum mechanical level. Given 
that a single realization of a BEC had a well-defined num- 
ber N, initially with all atoms in the ground state component 
|V)v(? = 0)) = \0,n g = N), then there is obviously no need 
for an extra number filtering device, as it is sharply defined by 
perfect preparation. However, the experimental reality usually 
is plagued with technical imperfections, day-to-day variations 
or finite temperatures ensembles. Thus, one should include 
the possible number uncertainty in a statistical description and 
use the density matrix 



N=0 



*n \y N {t))(VN(t)\, 



(18) 



where the distribution & accounts for the lack of informa- 
tion H3I1 . As usual, we obtain averages of observable with 
respect to the mixed state by a trace (...) = Tr{. . .p(f)} over 
all the Fock space. 
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FIG. 5: Evolution of the particle number distribution |i//JJ,(2r)| 2 
from an initially pure Fock state (solid vertical line) to a broadened 
distribution, exiting in the interferometer port e. In subplot a), we 
used |0,./V= 1720), which leads to a large number dispersion. In 
subplot b) we used the particle number |0,yV = 940) that defines the 
optimal working point of the interferometer in Fig. [4] 



The particle output in the e-channel of the interferometric 
number filter and its uncertainty are measured by averages, 
variances and volatility 

(n e {t)) = £ ^ N {y N (t)\n e \y N (t)) = 



N=0 



L^vE(a^«)i^(oi 2 , 

N=0 n=0 

nut °e{2T) 



(n e (2T)) 



(19) 



(20) 



We have now evaluated the time-dependent difference 
Schrodinger equation Eq. (TT3T > for the p — % — p interferom- 
eter for a range of initial particle numbers < N < 2500 in 
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the BEC. In particular, we have even relaxed the Raman-Nath 
approximation for the p — n — p beam-splitter sequence and 
used real rectangular pulses with duration T = 1.31 ms< T 
and total pulse areas Q.T such that p = p = 0.9, as before. All 
other parameters were as in Sec. IHI B~2l The results of this 
quantum mechanical simulation, in particular the mean out- 
put number (n e (2T)) and normalized number uncertainty s° ut 
are depicted in Fig. [4] They agree well with the homogeneous 
and inhomogeneous mean-field calculations. 

In Fig. [5] we show two special wavefunctions: one that cor- 
responds to the optimal input number N = 940 and one with 
N = 1720 that exhibits a large number dispersion after passing 
through the interferometer. The Gaussian nature of the final 
wavefunction can be explained from a semiclassical analysis 
of the difference equation [56]. Both wave functions can be 
identified also clearly as extrema of the response function in 
Fig. [4] From these two extreme examples of the evolution 
of pure initial Fock states it is obvious that a further con- 
volution with a non-ideally prepared ensemble, for example 
&n ~ ex p[~ (N — N) 2 /<j], will exhibit a similar dispersion 
response: a suppression of number fluctuation in the output 
channel if the optimal working point of the filter coincides 
with the most likely particle number of the input ensemble. 



IV. CONCLUSION 



In conclusion, we have analyzed the performance of a non- 
linear number filter for matter waves. This addresses the 
problem of technical shot-to-shot variation of particle number 
in current Bose-Einstein condensate experiments. A highly 
asymmetric p — it — p beam-splitter sequence is required to 
achieve optimal filtering performance. In the case of symmet- 
ric splitting or the absence of nonlinearity, no filtering is seen 
at all. This method is in direct analogy to a nonlinear fiber op- 
tics setup studied in if^L [tL FTotl . In detail, we have analyzed the 
situation for an homogeneous system and an inhomogeneous 
trapped gas within the Gross-Pitaevskii mean-field theory, as 
well as a quantum mechanical Josephson model, which ad- 
dresses complementary aspects and agrees well otherwise. 
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